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l. INTRODUCTION 


Ais BACKGROUND 

In the past most of the interest in ocean acoustics focused on open 
ocean acoustics where the ocean depth can be considered infinite, and 
hence interactions with the bottom not significant. However, there has 
been a recent shift in interest to shallow water acoustics which has in 
turn increased the interest in bottom interactions since they are so common 
in shallow water. Like the sea surface, the bottom is an effective reflector 
and scatterer of sound and acts to redistribute in the ocean above it a 
portion of the sound incident upon it (Urick, 1983, p.271). 

The ocean bottom can be modeled as a fluid-filled porous solid and 
hence general research into the acoustic properties of porous media is of 
interest. Specifically, the solid-fluid interaction in a fluid-filled porous 
solid can provide insight into bottom interactions. 

One way to characterize the microscopic solid-fluid interaction in a 
fluid-filled porous solid is through the use of the concept of complex 
effective mass density. It can be related to other familiar terms used in 
porous media acoustics, such as complex, or dynamic, tortuosity. In 
addition, other terms important in porous media research can be obtained 
from the complex effective mass density, such as effective fluid mass 


density and flow resistance. 


B. OBJECTIVES 

The objective of this research is to experimentally determine the 
complex effective mass density of air in a rigid porous solid over the 
frequency band of 10 to 1000 Hz. The experimental method being 
investigated involves using acoustic pressure and _e acceleration 
measurements to obtain the desired results. The theory is based on 


Newton's second law for a viscous fluid contained in a rigid porous solid. 


li. THEORY 


A. COMPLEX EFFECTIVE MASS DENSITY 
Newton's second law for a fluid, ignoring all sources of dissipation 


except shear viscosity, is: 


au = a i (2.1) 
a + p(u-V)u = -Vp + nVxVxu, 


where p is the fluid mass density, & is the fluid particle velocity, p is 
the fluid pressure, and is the fluid shear viscostty. For small 
displacements, such as occur during the propagation of a low-amplitude 
sound wave, the quadratic terms in Equation (2.1) can be neglected, 


resulting in the linear differential equation 


=p 


Ou - 
— = -Vp + nVxVxu. 
P pt+yn x 


(2.2) 


The corresponding linearized equation for a fluid-saturated rigid and 


stationary porous solid in terms of the average fluid velocity, (u) , is 


(Baker, 1986, pp. 334-335): 


Pay A) = - V1) - Rroeld) (2.3) 


where (p) is the average fluid pressure, Puy is the effective fluid mass 


density and Riew is the flow resistance. P,,y and Ry, are 


operationally defined by: 
\ 2 

5 Pag) 

2 

Ray) 


In addition, Py is always equal to or greater than the bulk fluid mass 


kinetic energy per untt fluid volume, 


rate of dissipation of fluid kinetic energy density. 


denstty, owing to the tortuostty of the pore geometry and to the 
nonuntformity of the microscopic velocity distribution across a pore. 


Equation (2.3) can be written as: 


au) 


eet Ryoyit) = - Vip), a 


or, assuming a time harmonic solution of the form elie », as: 


Au) _ (2.5) 
= Vip) . 


bets 





Equation (2.5) can now be written in terms of the average fluid 


acceleration, (a) : 
Bag (2) = - Vip) (2.6) 
where, 
Anew (2.1) 


Poy ™ Peg * jo. 


is the complex effective mass density, the quantity of interest. 


It is assumed that all quantities being measured have been averaged 
over the microscopic pore geometry of the porous solid, and since this 
research is only concerned with macroscopic phystal quantities, the 


brackets indicating such quantities will henceforth be omitted. 


B. EXPERIMENTAL METHOD 


An experiment was designed and conducted in which a quantity of air 
is made to oscillate back and forth as a plug through a porous solid 
sample such that the pressure is zero at the midplane of the sample and 
equal in magnitude but opposite in phase at the ends. Figure 1 shows the 


setup used for the experiment. 


Ze 


- 


ne 





Figure 1 Experimental setup. 


The complex mass density of the air can be obtained from pressure 


and acceleration measurements taken at the ends of the sample as follows. 


If it assumed that the only motion is one-dimensional motion along the axis 


of the sample, Equation (2.6) can be rewritten as: 


= _Vp (2.8) 


Because the air is oscillated at very low frequencies, the gradient in 
pressure along the sample can be approximated as the ratio of the 


pressure difference along the sample to the length of the sample, 


Vp « a (2.9) 


where Pp, and p, are the pressures at the right and left ends of the 


sample, respectively. 

Equation (2.9) assumes that the wavelength of sound in air, A, is much 
much greater than L. If this is not the case then a finite wavelength 
correction can be applied to Equation (2.9). Although, A>>L_ at all 
frequencies in the experiments conducted, a finite wavelength correction 
was applied to improve the accuracy of the results. The finite wavelength 
correction factor for the empty tube experiments is discussed in Chapter 
IV and for the porous solid sample experiments in Chapter V. 

Similarly, the acceleration of the air within the sample can be 


approximated as the average of the acceleration at the ends of the sample, 


a,+a 
‘: iam (2.10) 


where @, and @, are the acceleration at the right and left ends of the 


sample, respectively. 
Hence, using pressure and acceleration measurements taken at the 
ends of the sample the complex effective mass density of the air can be 


approximated by 


By * -= 2. (2.11) 


lil. DESCRIPTION OF EXPERIMENT 


A. APPARATUS 


The apparatus shown in Figure 2 was built to hold an air-filled 
porous solid sample. The apparatus holds the porous solid sample within 
a cylindrical sample tube which is capped on each end by a speaker end 


cap assembly (speaker end cap). 





Figure 2 Apparatus used in complex mass density experiments. 


ie Speaker End Caps 
The speaker end caps, shown in Figure 3, consist of two identical 
moving-coil loudspeakers, each of which are attached to an acrylic 
mounting ring. The loudspeakers are four-inch diameter woofers 
manufactured by Radio Shack, Model 40-1022B (Tandy Corporation, Fort 
Worth, TX.). Modifications were made to each speaker end cap in order to 
take pressure and acceleration measurements. 

To measure the pressure at each end of the sample, microphones 
were attached to the speaker end caps. The microphones are Radio Shack 
electret condenser microphones, Model 270-090, which were attached to two 
stainless steel wires posttioned over the cone. In addition to supporting 
the microphone, the stainless steel wires also served as the electrical 
connections. 

An accelerometer was attached to each speaker cone to measure 
acceleration at each end of the sample. An Endevco ISOTRON™ piezoelectric 
accelerometer, Model 25R, (Endevco Corporation, San Juan Capistrano, CA.) 
was selected for its extremely small size (approx 2 mm square) and weight 
(140 mg). The small dots on the surface of the each speaker cone were 


used in the accelerometer calibration procedure described below. 


2. Sample Tube 
A sample tube made of cast acrylic was chosen to hold the rigid 
porous solid sample. A sample tube 4.0 inches in length with an inner 
diameter of 3.5 inches and an outer diameter of 4.0 inches was selected. 


The length of the sample tube was chosen to minimize finite wavelength 


effects throughout the entire frequency range of interest. The inner 
diameter of the tube was chosen to closely match the effective radiating 


diameter of the loudspeakers used in the speaker end caps. 





Figure 3 Speaker end cap assemblies. 


3. Sensor Cahbration Procedures 
The accelerometers were calibrated by the manufacturer and a 
rated sensitivity was provided with each accelerometer. However, tt was 
necessary to verify that the accelerometers were, in fact, measuring the 
acceleration of the air slug passing through the sample, and not just the 


acceleration at only one point on the speaker cone. 
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The condenser microphones selected as the “pressure sensors” 
were not calibrated by the manufacturer. Therefore tt was necessary to 


calibrate them. 


a. Accelerometer Calibration 

As mentioned above, the accelerometers were calibrated by 
the manufacturer; therefore, a complete calibration was not necessary. It 
was, however, necessary to verify that the acceleration being measured was 
the same as the acceleration at any other point on the speaker cone. The 
idea was to compare the output of the accelerometers to another 
independent measure of acceleration. 

It was not desirable to attach another device to the speaker 
end caps to provide the second independent measure of acceleration since 
the additional weight of the second device might affect the results. 
Instead, another option was chosen, which was to measure the velocity of 
the cone by a noncontact means using a laser doppler velocimeter and to 
compare the results to the acceleration as measured by the accelerometer. 

The relationship between acceleration and velocity is given 


by 


=@, (3.1) 


where: 
u = veloctty, 


acceleration. 
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In addition 


u- Ae and a = jwAe™, (3.2) 


therefore 


: = jo = janf. (3.3) 


Hence the ratio of acceleration to velocity is a linear function of frequency. 

As mentioned above, the velocity of the speaker cone was 
measured using a laser vibrometer, a device which measures velocity by 
comparing the doppler shift in a light beam which has been reflected off 
the object whose velocity is desired. Therefore, the ratio of accelerometer 
output to vibrometer output should be a linear function of frequency. 

Equation (3.3) above can be stated in terms of the outputs 
of the accelerometer and laser vibrometer as: 


{accelerometer (mV) Wocceleromeser Vig)" mis"/g)\ . asf, 
[vibrometer voltage (mV)jivelocity range (m/s{V)) 





or, 
[accelerometer voltage (mV)) _ [velocity range (mj{V)}leccelerometer_sensitivity( Vig) uf . 
|: toromen: vokage (mY) [gravity (m/s/g)) 


It can easily be seen this is an equation of a straight line when expressed 


in the form: 


[accelerometer voltage (mV)] _ S. 2nf. 
[vibrometer voltage (mV)] 9 


Hence, if the ratio of accelerometer output voltage to laser vibrometer 
output voltage is plotted as a function of frequency, the curve should be 


as shown in Figure 4. 
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Figure 4 Theoretical comparison of accelerometer and vibrometer output. 


One major benefit of using the laser vibrometer for 
measuring velocity is that it is capable of measuring the velocity at any 
point on the cone. Hence, if tt could be shown that the ratio of 
accelerometer output to laser vibrometer output was a linear function of 
frequency, independent of the location on the speaker cone, the 
verification would be accomplished. 

The setup shown in Figure 5 was used to verify that the 
accelerometers were, in fact, measuring the acceleration of the air slug 
passing through the sample, and not just the acceleration at only one point 


on the speaker cone. 
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Figure 5 Accelerometer - 


vibrometer comparison setup. 


Laser vibrometer output at several points on each speaker 
cone was compared to accelerometer output. 


Small dots of reflective 
material were placed on the surface of the speaker cone at each point to 
reflect the laser light beam. 
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The ratio of accelerometer output to laser vibrometer output 
for the right and left speaker end caps are shown in Figures 6 and 7, 
respectively. The theoretical lines were derived using a calibration factor 
of 25 mm/s/v for the laser doppler vibrometer and the _ individual 
accelerometer calibration constants supplied by the manufactuer. In both 
cases the ratio of accelerometer output to vibrometer output is essentially 
a linear function of frequency. The graphs shown are for a single 
comparison point on each speaker end cap; graphs for other points are 


similar. 
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Figure 6 Ratto of accelerometer output to vibrometer output. 
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Figure T Ratio of accelerometer output to vibrometer output. 


Since the ratio of accelerometer output to vibrometer output 
was found to be a linear function of frequency, independent of the location 
on the speaker cone, it is concluded that the output of the accelerometers 
is, in fact, a measure of the acceleration of the air slug passing through 
the sample. 

b. Pressure Sensor Calibration 

As previously mentioned, the condenser microphones selected 
as the "pressure sensors" were not calibrated by the manufacturer; 
therefore it was necessary to calibrate them prior to use. Two independent 
calibration methods were chosen, first, a comparison calibration, and 


second, with a pistonphone. 
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(1) Comparison Calibration. The condenser microphones 
were first calibrated using a comparison method. In this method, a 
comparison is made of the output of the condenser microphone to the 
output of a calibrated microphone while both are being subjected to the 
same sound field. Then, given the known sensitivity of the calibrated 
microphone in dB re lv / Pa, the output of the calibrated microphone in 
dB re iv / Pa, and the output of the condenser microphone in dB re iv 
/ Pa the sensitivity of the condenser microphone could be determined as 


follows: 


unknown sensitivity level = known sensitivity level + difference in output levels (3.4) 


The setup shown in Figure 8 was used to accomplish 
the comparison. The moving-coil loudspeakers mounted in the speaker end 
caps were used to produce the sound field. A 1/8 inch diameter Bruel & 
Kjaer condenser microphone, Type 4138 (Bruel & Kjaer, Foster City, Ca.), 
with a sensitivity of - 68 dB reilv / Pa (0.40 mV / Pa) was the calibrated 
microphone used for the comparison. A Hewlett-Packard Impedance 
Analyzer, Model 4194A (Hewlett-Packard Corporation, San Jose, Ca.), was 
used to measure and record the ratio of the output of the “pressure 
sensor” to the output of the B&K microphone, since it could do so while 
sweeping the drive signal to the loudspeakers over the desired frequency 


band, 10 to 1000 Hz. 
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Figure 8 Comparison calibration setup. 


As can be seen in Figures 9 and 10, the sensitivity of 
both "pressure sensors" was a function of frequency. A least squares fit 
of the data was used to determine this relationship with respect to 


frequency. 
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Figure 9 


“pressure sensor” / B&K Microphone (d8) 


Figure 10 


right end cop pressure seneor 


—~-~-— least equores fit 24.13 — 0.0015 @ freq 
ratleo(250 Hz) = 23.7 dB 





0 100 200 300 400 500 600 700 800 900° 1000 
frequency (Hz) 
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left end cap pressure sensor 


-~--— least squores fit 25.55 — 0.00035 @ freq 


ratio(250 Hz) = 23.4 dB 
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Ratio of "pressure sensor” output to B&K microphone output. 
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The resulting ratios of the output of the "pressure 


sensors’ to the output of the B&K microphone were: 


right pressure sensor / B&K = 24.1 - 0.0015 f dB, 


left pressure sensor / B&K = 23.5 - 0.00035 f dB, 


where f is frequency in Hz. 


Given the above ratios and the sensitivity of the B&K microphone the 


unknown sensitivities were determined using Equation (3.4): 


right pressure sensor sensitivity = - 43.9 - 0.0015 f dBre 1V / Pa, 


left pressure sensor sensitivity = - 44.7 - 0.00035 f dBre 1V / Pa. 


(2) Calibration by Pistonphone. The second method used to 
calibrate the ‘pressure sensors" was by pistonphone. A B&K Pistonphone, 
Model 4220, with a calibrated output of 123.9 dB re 20 uPa at 250 Hz was 
used. Had the frequency response of the "pressure sensors’ been flat 
this method could have been used to verify the comparison calibration over 
the entire frequency band of interest. Unfortunately this was not the 
case; therefore, the pistonphone could only be used to verify the 
comparison calibration at 250 Hz. 

The setup shown in Figure 11 was used. An adapter 
was needed since the pistonphone was designed to calibrate B&K 
microphones, which have a slightly smaller case than the Radio Shack 


condenser microphones being used. 
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Cahbration by Pistonphone. 


Figure 11 


The sensttivtties found using the B&K Pistonphone were: 


- 43.7 dB re 1V / Pa, 


right pressure sensor sensttivity 
- 44.6 dB re 1V / Pa. 


left pressure sensor senstivity 
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At 250 Hz the sensitivities of the ‘pressure sensors" as found by the 


comparison calibration were: 


right pressure sensor sensitivity = - 44.3 dB re 1V / Pa, 


left pressure sensor sensitivity = - 44.7 dB re 1V / Pa. 


These values are within 10% (1dB) of the sensitivities found using the B&Kk 
Pistonphone which would indicate the comparison calibration is valid, at 


least at 250 Hz. 


B. EQUIPMENT SETUP 

The equipment setup shown in Figure 12 was used for the experiments 
described below. The basic idea was to drive the moving-coil loudspeakers 
so that a null in pressure would be obtained at the center of the sample 
tube. Measurements of pressure and acceleration needed to be made at a 
set of frequencies chosen within the desired band of 10 to 1000 Hz. 

In order to obtain a null at the center, tt is necessary that the 
loudspeakers be driven out of phase. To prevent the phase from 
wandering, a method of locking the two phases together had to be found. 
The solution was to use two HP Function Generators, Model 3314, which 
were connected such that one was acting as a master and the second as 
a slave. The sync output of the master was connected to the trigger input 
of the slave. The major limitation to this setup was that it could only be 
used at 50 Hz and above due to limitations of the generators. Another 


minor limitation was that, because of their 50 ohm source impedance, the 
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HP 3314's could not provide a stable enough output voltage when loaded 


down by the loudspeakers. Therefore, the output of each HP 3314 was 


passed through an audio frequency amplifier, a tTechron Model 7520 Power 


Supply Amplifier (Crown international, Inc., Elkhart, In.). 
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Figure 12 Equipment setup. 
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The output of each HP 3314 was monitored using a Kikusui Dual 
Channel Oscilloscope, Model DSS5020. This setup allowed the relationship 
between the phase of the two loudspeaker signals to be monitored on the 
same display. 

Measurements were taken using a Stanford Research Systems Lock-in 
Amplifier, Model SR530 (Stanford Research Systems, Sunnyvale, Ca.). A 
lock-in amplifier was chosen since both magnitude and phase measurements 
needed to be taken. A second Kikusui Oscilloscope was used to monitor the 
signal being sent to the lock-in amplifier. 

A connection box was designed and built to allow selection between 
each of the sensors. The output of the connection box was connected to 
the input of the lock-in amplifier and the oscilloscope. 

All trigger inputs were connected to the sync output of the master 


HP 3314 for a single reference signal. 
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iV. EMPTY TUBE EXPERIMENTS 


A. INITIAL EMPTY TUBE EXPERIMENT -~- THEORY 
For the first experiment, an empty sample tube was used, as shown 


in Figure 13. 











Figure 13 Empty tube experimental setup. 


The objective of this experiment was to test the apparatus and equipment 
setup on a known quantity, namely air in an empty sample tube. The 
density of air at 20°C is 1.21 kg/m* (KFCS, 1982, Appendix 10). 

The complex effective mass density of the air contained in the empty 
sample tube, as defined in Chapter !!, is the negative of the ratio of the 
gradient in pressure along the sample tube to the acceleration of the air 


within the sample tube, that is 


im =\p (4.1) 


Zo 


As previously discussed above, the gradient in pressure along the 
tube can be approximated as the ratio of the pressure difference along the 


sample tube to tts length, 


Pre-~ Pp _ = , 
Vp a LR FL, Vp A (4.2) 
L 
where Pp, and p, are the pressure at the right and left ends of the 
sample tube, respectively. 


The acceleration of the air within the sample tube can be approximated 


as the average of the acceleration at each end of the sample tube, 


+ ms 
a Rot aa, (4.3) 


where @, and @, are the acceleration at the right and left ends of the 


sample tube, respectively. Hence, Equation (4.1) can be written as 


2P-m, Wp 


‘ 4.4 


Poy 

The calculation of Poy using Equations (4.2) through (4.4) appears to 

be rather simple; however, corrections need to be made to the above 
equations for the effects of sound propagation within the air column, and 
the slight difference in the cross-sectional area of the sample tube and the 
loudspeakers at either end. These corrections are discussed in the 


following subsections. 
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1. Finite Wavelength Correction Factor 

Equation 4.2 assumes the air in the tube moves as a solid plug 
and as a result, that the pressure drop across the sample tube is 
distributed in a linear fashion as shown in Figure 14a. Unfortunately, 
plug flow is the limiting case when the wavelength of the oscillation of the 
air in the tube is infinite. For the cases of interest the wavelength is 
finite and the pressure drop is distributed approximately as shown in 
Figure 14b. (The pressure distribution shown in (b) is for a pressure 


wave with a wavelength of 2L and at one specific point in time.) 





Figure 14 Pressure distribution within sample _ tube. (1) infinite 
wavelength, (2) finite wavelength. 


The acoustic pressure and volume velocity at any point within 
the sample tube can be related using the concept of a velocity potential 
(¢). The following relationships exist among acoustic pressure (p), particle 
velocity (u), and velocity potential (¢) for plane wave propagation in one 
dimension (generalized to the case of porous media from Kinsler,et a/., 


1982, pp. 104-107): 


p=-j@OPpy®, (4.5) 
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and 


ur, (4.6) 
For the fintte wavelength case the pressure distribution is 


p(x) = Do sin kx , (4.7) 


where & is the wave number of the pressure wave, and hence 


— - 2p, sin KL/2 
oe 


. (4.8) 


Substituting Equation (4.7) into Equation (4.5) and solving for ¢ yields 





sin kx 
@(x) = - Po (4.9) 
Jj @ Poy 
Equation (4.6) then becomes 
K Dp, COB kx 
syn! eee (4.10) 
J @ Poy 
Therefore, the particle acceleration is 
Po 
a(x) =-k cos kx , (4.11) 
off 
and the acceleration of the air at each end of the sample tube is 
Po 
a(+L{/2) = - k — cos kLf/2. (4.12) 
of 
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Therefore, 


Bo a ec kL/2 . (4.13) 
eff 


Equations (4.8) and (4.13) can be combined to yield 


_ Sp. ae 
= Par| zp |" (4.14) 


re WE : ae (4.15) 
a 


where 


(4.16) 


Fes = |e Ee) 


kL/2 


is the finite wavelength correction factor. 


2. Acceleration Correction Factor 
As previously stated, the acceleration of the air within the sample 
tube is taken to be the average of the acceleration at each end of the 
sample tube. However, the acceleration being measured is the acceleration 
of the speaker cones. Unfortunately, the acceleration of the speaker cones 
and the acceleration at the ends of the sample tube are not equal because 
the cross-sectional area of the sample tube and loudspeaker are not quite 


equal. However, the mass flow rate of the air at the surface of the 
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speaker cone is equal to the mass flow rate of air at the end of the sample 


tube. Hence, 


Sione Acont™ Frube Aeube (4.17) 
where 
Gone = acceleration of the surface of the speaker cone, 
Giube = acceleration of the air at the end of the sample tube, 
A..ng = Cross-sectional area of the speaker cone, 
A.4, = cross-sectional area of the end of the sample tube, 


Solving Equation (4.17) for the desired quantity, @.,. , yields: 








A Dalia 
Cone Cone 
Oinpe * Scone | ~ Scone | (4.18) 
A Disre 
where 
Done = Aiameter of the speaker cone (8.3 cm.), 
Dab = diameter of the sample tube (8.9 cm.). 
Hence, Equation (4.3) can be written as: 
Q@a2%+@ a, +@ 
jg lh (4.19) 


2 2 


where ay and a, are the acceleration of the right and left speaker 


cones, respectively, and 


D 2 
acf = aa = 0.870 (4.20) 
Deabe 
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is the acceleration correction factor. Equation (4.15) then becomes 


ee el 
off Gy +O, ac fwof (4.21) 
2 


B. INITIAL EMPTY TUBE EXPERIMENT - RESULTS 

Measurements were taken over the frequency range of 50 to 500 Hz 
in 10 Hz steps and over the frequency range of 500 to 1000 Hz in 50 Hz 
steps. At each frequency the amplitude and phase of the slave HP 3314 
was adjusted to obtain a pressure null at the center of the sample tube. 
it was assumed when this condition was obtained that the two loudspeakers 
were moving exactly out of phase. Once the null was obtained, pressure 
and acceleration measurements were taken at each end using the SR Lock- 
in Amplifier. Both amplitude and phase measurements were taken since a 
complex quantity was being calculated. 

A computer program was written to compute the complex effective mass 
density at each frequency using Equation (4.21). Then, since 


Puy i Puy =p nee 9 (4.22) 


the effective fluid mass density is 
Poy * Re[p,,] 9 (4.23) 
and the flow resistance is 


Rasy = ~ Olm[P yl - (4.24) 


A copy of the program is attached as Appendix A. 
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Although the results, as shown in Figures 15 and 16, were initially 
encouraging, the steady increase in both the density and flow resistance 


at frequencies above 700 Hz was unexpected and thus warranted further 


investigation. 


denalty = 1.21 kg / ma3 
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Figure 15 Initial empty tube results - effective fluid mass denstty. 
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Figure 16 Initial empty tube results - flow resistance. 
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Careful examination of the data suggested a possible problem with the 
pressure phase measurements. Although the condenser microphones used 
as "pressure sensors’ were calibrated with respect to pressure magnitude 
there was no calibration made with respect to pressure phase. A 
comparison with respect to pressure phase was made with the calibrated 
B&K microphone used to calibrate the sensors, since itt has a flat phase 
response. The results, shown in Figure 17, indicate that this is a potential 


source of the problem. 


phese angie relative to BAK 


retative phase angie (deg) 
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Figure 17 Pressure sensor comparison to B&K with respect to phase. 


Attempts were made to correct for this problem; however, itt was 
decided that since more accurate sensors were available, modifying the 
apparatus would be a simpler solution. In addition, a revised experimental 
procedure was developed to correct for another potential source of error 


due to air compression. 
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C. MODIFICATIONS TO APPARATUS 

The left end cap pressure sensor was replaced with a PCB 
Piezotronics microphone, Model 103A (PCB Piezotronics, Inc., New York, NY). 
This microphone is rated to 25 psi with a sensitivity of 2080 mV/psi or 
0.3017 mV/Pa and was supplied with a calibration curve good from 0 - 3 
psi. The modification to the left speaker end cap changed the effective 
length of the sample to 92.4 mm vice 102 mm. 

The center pressure sensor was also replaced with a PCB 
Piezontronics microphone, Model 103A, with a rated sensitivity of 1913 
mV/psi or 0.2775 mV/Pa. The right end cap pressure sensor was not 
replaced for reasons which will become clear when the revised experimental 


procedure is discussed. 


D. REVISED EMPTY TUBE EXPERIMENT - THEORY 

The acceleration correction factor previously developed only corrected 
for the difference in the cross-sectional areas of speaker cone and the 
tube. There was however, another potential source of error that had not 
been considered, and that was the acceleration loss between the surface 
of the speaker cone and the end of the sample due to the compressibility 
of air. Hence, either an additional correction factor needed to be applied 
to Equation (4.21) to correct for that loss or the procedure needed to be 
revised to eliminate the loss. It was decided that the best approach was 
to revise the procedure to eliminate the loss. 

Instead of adjusting the slave HP 3314 to obtain a null in pressure 


at the center of the sample tube, the slave HP 3314 was adjusted to obtain 


a null at the right end of the sample tube. The resulting pressure drop 
across the sample is distributed as shown in Figure 18. The wavelength 
of the pressure wave is finite in all cases. (The pressure distribution 
shown is for a pressure wave with a wavelength of 4L and at one specific 


point in time. ) 
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ieee ne 


Figure 18 Pressure distribution within sample tube with pressure null at 
the right end of the sample. 


As aresult, acceleration loss between the surface of the speaker cone 
and the end of the sample due to the compressibility of air is identically 
zero (since the acoustic pressure there is zero) and the acceleration of the 
air within the sample tube is the acceleration at the right end of the 


sample tube. Therefore, 
a= @,7 a - acf . (4.25) 


The gradient in pressure along the sample tube can _ still be 
approximated as the ratio of the pressure difference along the sample tube 
to its length; however, since the pressure is zero at the right end of the 


sample the pressure difference along the sample is just the pressure at 


3° 


the left end of the sample. Therefore, Equation (4.2) can be written as 


~P, 


Vp = = Vp. (4.26) 





In addition to the changes listed above, the finite wavelength 
correction factor as defined by Equation (4.16) needed to be changed since 
the pressure drop across the sample tube was no longer distributed as 


shown in Figure 14. 


1. Revised Finite Wavelength Correction Factor 
As previously stated, the acoustic pressure and volume velocity 
at any point within the sample tube can be related using the following 


relationships for plane wave propagation in one dimension: 


P= -JXpy®, (4.27) 
and 
d@® 
— ° 4.28 
etl Ce 


For the case the shown in Figure 18, the pressure distribution is 
P(x) = - Pp Sim kx, (4.29) 


where k is the wave number of the pressure wave, and hence 


= _ Po Sin KL 


Vp = (4.30) 
P L 
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Substituting Equation (4.29) into Equation (4.27) and solving for ¢ yields 


sin Ax 
@(x) = noe : (4.31) 
jo Puy 
Equation (4.28) then becomes 
cos kx 
u(x) = Po (4.32) 
1 Dy 
Therefore, the particle acceleration is 
Po 
a(x) =k cos kx , (4.33) 





eff 


and the acceleration of the air at the right end of the sample tube is 


P 








a, = a(0) =k —. (4.34) 
off 
Therefore, 
a, -acf=k Po (4.35) 
off 
Equations (4.30) and (4.35) can be combined to yield 
_ Up sin kL 
a, - acf Pup kL |. Se 


Solving Equation (4.36) for the complex effective mass density then gives 


= = __Vp e _i 
Pag a, ‘acf fwef P 
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where 


(4.38) 


‘tied oye H 


is the revised finite wavelength correction factor. 


E. REVISED EMPTY TUBE EXPERIMENT - RESULTS 

Measurements were again taken over the frequency range of 50 to 500 
Hz in 10 Hz steps and over the frequency range of 500 to 1000 Hz in 50 
Hz steps. However, this time the amplitude and phase of the slave HP 3314 
was adjusted at each step to obtain a null in pressure at the right end of 
the sample tube. Once the null was obtained, pressure and acceleration 
measurements (amplitude and phase) were taken at the left and right ends 
of the sample tube, respectively, using the SR Lock-in Amplifier. 

A computer program was written to compute the complex effective mass 
density at each frequency using Equation (4.37). Then, since 


-j ee , (4.22) 


Pay * Pet 


the effective fluid mass denstty is 


Pag = RelP,yl » (4.23) 


and the flow resistance is 
Ry = ~ olm[A,,! . (4.24) 


A copy of the program is attached as Appendix B. 
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The revised empty tube results for effective fluid mass density and 
flow resistance are shown in Figures 19 and 20 respectively. The change 
in procedure had a dramatic effect on the effective fluid mass density 
above 700 Hz. It also had an effect on the flow resistance, but not so 


dramatic. 
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Figure 19 Revised empty tube results - effective fluid mass denstty. 





0 100 200 300 406 8 800 @00 700 eno 600 0=—_« 1000 
frequency (Hz) 


Figure 20 Revised empty tube results - flow resistance. 
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it is apparent from the results oe in Figure 19 that the objective 
of this experiment was achieved, namely, to test the apparatus and 
equipment setup on a known quantity, i.e. air in an empty sample tube. 

At this point the interpretation of the results for flow resistance 
shown in Figure 20 are not of a major concern, since flow resistance is not 
really a relevant quantity for air in an empty tube. In addition, a 
thorough discussion of flow resistance requires the introduction of 
concepts, such as viscous penetration depth, which are much more relevant 
to air in a porous solid. Therefore, these concepts will be introduced in 
the next chapter and a complete discussion of flow resistance will follow. 
The results of both empty tube and porous solid experiments will be 


included in that discussion. 
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V. POROUS SOLID SAMPLE EXPERIMENTS 


A. POROUS SOLID SAMPLE SELECTION 
The measure of the influence of fluid shear viscosity on a fluid 
saturated porous solid subject to elastic wave propagation is the viscous 


penetration depth, defined as: 


% 
5 isc 23 ’ (5.1) 


op 





where p is the fluid mass density, is the fluid shear viscosity, and W is 
the frequency. The ratio of viscous penetration depth to a characteristic 
pore size is one of the most fundamental dimensionless parameters that 
determines the nature of elastic wave propagation in a fluid saturated 
porous solid (Baker,1986, p. 27). 

if the viscous penetration depth is much less than a characteristic 
pore size, then the motion of the fluid is very similar to that of an 
inviscid fluid, and independent (of the solid) sound propagation within the 
fluid is possible. in the porous solid literature this case is referred to as 
the high frequency or boundary layer flow limit. 

if the viscous penetration depth is much greater than a characteristic 
pore size, then the fluid in the pores is essentially frozen to the solid and 
sound cannot propagate within the fluid. This case is referred to as the 
low frequency or Poiseuille flow limit. 

Both cases are of interest and were investigated. The measure of 


characteristic pore size used was mean pore diameter since this was the 
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only measure the porous solid sample manufacturer (General Polymeric, West 


Reading, PA.) could work with. 


Three samples with mean pore diameters 


of 150, 300, and 500 microns were selected so that over the frequency band 


of interest, 50 to 1000 Hz, both types of flow could be observed and the 


transition between the two types of flow would occur at several different 


frequencies. 


Figure 


20 shows the estimated transttion from Poiseuille flow to 


boundary layer flow (plug flow) for a 300 micron sample. 


viscous pentration depth / pore rodius 


mean pore diameter 300 microns 


onsen 
ras. 


o ea 


7 ° 
a * ‘* 


ae 





10! 
‘Sey ttt ee Seledae pe wigs Ve “sey 
tere at ere 70 ee ae set 
: i i 1 holo ts alae ae : tom 
Secretar ee pier esses root oOo ae ele 

; ae ic a ae ; ee i 
Poisevitle flow : 
10° : 
Seat ce ter edece sce deerectes apiece eee 46 
peo bee sce 
s & & 
107! ° s . e e s * . ° 
101 10? 


frequency (Hz) 


105 


Figure 20 Estimated flow transition for a 300 micron sample. 


B. THEORY 


The procedure developed for the revised empty tube experiment was 


used since it avoids corrections for the acceleration loss due to the 


compressibility of the air between the speakers and the ends of the 
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sample. The empty sample tube was, however, replaced with a sample tube 


containing a porous solid sample as shown in Figure 21. 
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Figure 21 Porous solid sample experimental setup. 


The slave HP 3314 was again adjusted to obtain a null in pressure at 
right end of the sample tube. Therefore, the complex effective mass 


denstty of the air contained in the sample can be computed using 


= Pp @ 1 a 
Pat GQ, -acf fwef C2) 
where 

—— Pr 
V = —_ a 9 (5.3) 

P L 

sin kL 

= —— 5.4 
fwef came (5.4) 


and Gy is the acceleration of the right speaker cone. The acceleration 


correction factor, acf , used for the empty tube experiments; however, 


needs to be modified. 
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1. Modified Acceleration ete Factor 
The method used to develop the acceleration correction factor for 
the empty tube experiments can again be used; however, this time the mass 
flow rate of the air at the surface of the speaker cone and the mass flow 


rate of air at the end of the sample tube are related as follows: 


Scone eons = Gsemple asm pie P. sample (5.5) 
where 
Gone = acceleration of the surface of the speaker cone, 
Gsample = acceleration of the air at the end of the sample, 
A.ong = Cross-sectional area of the speaker cone, 
Aseusie = cross-sectional area of the end of the sample, 
eamnte = porosity of the sample. 


The additional parameter, P , the porosity of the sample, is defined 


sample 
as the ratio of the volume occupied by the fluid to the total sample volume. 
It is required since the cross-sectional area of the fluid is reduced by the 


frame of the porous material. 


Solving Equation (5.5) for the desired quantity, Gample yields: 








2 
Byaaotg”™ Fcons 2 a = Gone * | i oo | (5.6) 
sample sample 
where 
D. une = diameter of the speaker cone (8.3 cm.), 
D, ample = diameter of the sample (8.9 cm.). 
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Hence, 


D sample 


D a 
acf = ie _ - + 0870 (5.7) 
P 
is the modified acceleration correction factor. The acceleration correction 
factor used for the empty tube experiments, Equation (4.20), can be 


obtained from Equation (5.7) by setting the porosity equal to one. 


a. Porosity Measurements 
The porous solid samples were constructed of glass spheres 
with a density, 9, of 2.4 g/om* bonded together to form a slug that was 
pressed inside of an empty sample tube. The porosity of each sample was 
determined by comparing the expected mass of a solid glass sample of 
identical size, constructed with glass of the same denstty, and the actual 
mass of the porous sample. 


The expected mass of a solid glass sample is 


po -volume =p-° 57h (5.8) 


where 
d= diameter of the sample (8.89 cm), 


height (length) of the sample (10.38 cm). 


Hence, the expected mass of a solid glass sample is 1546 g. 
The porosity of a given sample is a function of the ratio of 
the actual mass of the porous solid sample to the expected mass of a solid 


glass sample, that is 
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p «1 - —2ctual mass 


(5.9) 
expected mass 


The results of the porosity measurements are summarized in Table 5.1. 


Table 5.1 POROSITY MEASUREMENTS 


microns) x) 
pa | 50 | t002.3 | ons 
2 | 00s | ce | 
3s ft s00 873.5 | 
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2. Estimating kL in a Rigid Porous Solid. 

For the empty tube experiments, the wave number is the same as 

for free space, k = & with c * 340 m/s. However, for these experiments 
the tube was filled with a rigid porous solid, and the speed of the 
pressure wave was not known a priori. Therefore, a method had to be 
developed to determine the wave number. 
The method chosen was to calculate the wave number based on 

pressure measurements taken at the left end and the center of the sample 
tube as follows: 


The pressure at the left end of the sample is 


P, = P(-L) = ~p, sin k(-L) = py sin kL , : (5.10) 
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which can be written as 


P, = 2p, sin EB : cos ee] (5.11) 


and the pressure at the center of the sample, Dp, , is 


2 = (4) “Py sn K}-5] = py sin et 4. (5.12) 


Then taking the ratio of the two pressures, 


Pi 


= 2 cos = , (5.13) 
Pe 2 


Solving Equation (5.13) for kL in terms of the two pressures then gives 


kL = : cos”! 4 (5.14) 


This estimate was used in Equation (5.4) for the finite wavelength 


correction factor, fwef . 


Cc. RESULTS 

Measurements for each of the samples were made over the frequency 
range of 50 to 500 Hz in 10 Hz steps and over the frequency range of 500 
to 1000 Hz in 50 Hz steps. The amplitude and phase of the slave HP 3314 
was adjusted at each step to obtain a null in pressure at the right end of 


the sample tube. Once the null was obtained, pressure and acceleration 
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measurements (both amplitude and SRacey were taken using the SR Lock-in 
Amplifier, pressure at the center and the left end of the sample tube, 
acceleration at the right end of the sample tube (speaker cone). 

A computer program was written to compute the complex effective mass 
density of the air contained in the sample at each frequency using 


Equation (5.2). Then, since 


(5.15) 


j stem 
w 


bug = Poy = 


the effective fluid mass density of the air contained in the sample is 


Pay = Rel), (5.16) 


and the flow resistance is 
Raow = - wim[p). (5.17) 


A copy of the program is attached as Appendix C. 

Plots of the effective mass density and flow resistance of the air 
contained in the 150 micron sample are shown in Figures 22 and 23 
respectively. Similar plots for the 300 micron and 500 micron samples are 
shown in Figures 24 through 27. Included on each plot of the effective 
mass denstty for reference is the effective mass denstty of air contained 


in an empty sample tube at 20° C. 
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x — 150 micron sample data 
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Figure 22 Effective mass density - 150 micron sample. 
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Figure 23 Flow resistance - 150 micron sample. 
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x — 300 micron sample dato 
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Figure 24 Effective mass density - 300 micron sample. 
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Figure 25 Flow resistance - 300 micron samptke. 
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Figure 26 
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Effective mass density - 500 micron sample. 
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Figure 


frequency (Hz) 
27 Flow resistance - 500 micron sample. 
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D. INTERPRETATION OF RESULTS 
¢ 


1. Calculation of Puy using Biot’s Theory 


Using Biot's theory, the complex effective mass density can be 
expressed as follows (Baker, 1986, pp. 44, 48): 
Buy = &(W) py = ap, + Roe F(x) , (5.18) 
jo 
where @ is the (intrinsic) tortuosity, Py is the mass density of the bulk 
fluid, Ryo is the DC flow resistance, and F(x) is Biot's universal complex, 
frequency dependent correction function, whose purpose is to make the 
frequency dependence of buy come out right. F(x) is defined by (Baker, 
1986, p. 53) 


F(x) = 22 


F(x) = 4 Iq)’ (5.19) 


where 


3 R,- |-3 
z= J7(x) and x= 6 —s “. (5.20) 


& W py 





and 6 is the Biot "structural factor", a dimensionless parameter which 
depends only upon pore shape. The three parameters @ , Rp , and 
& completely characterize the microscopic pore geometry in Biot's theory. 
Although the tortuosity @ may be obtained experimentally, it is 
found to depend primarily on the porosity P in a nearly predictably 


fashion which for the case of a porous solid constructed of a random 
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collection of spheres is adequately represented by (Sen et. al. 1981, 


Johnson and Sen 1981) 


a= Pp-* . (5.21) 
Calculations of tortuosity for each of the samples using Equation (5.21) are 


summarized in Table 5.2. 


Table 5.2 TORTUOSITY CALCULATIONS 





The structural factor 6 depends weakly upon the cross-sectional 
shape of the pores. Baker suggests a value of 4 for 6 for a porous solid 
constructed of a random collection of spheres (Baker, 1986, pp. 291 -292). 

The DC flow resistance is proportional to the square of the ratio 
of the viscous penetration depth to a characteristic pore size (Baker, 1986, 


p. 45). Hence, 


q 


CO ——————— (5.22) 
(pore sizey 


Roc 


where Is the fluid shear viscostty. Therefore, the low frequency limit 
of the flow resistance should decrease with increasing pore diameter. This 


is clearly the case as shown in Figure 28. 


5S 


Rfiow (kg/s /m-3) 





© 100 200 300 400 800 600 #700 800 £900 1000 
frequency (Hz) 


Figure 28 Flow resistance - 150, 300, and 500 micron samples. 
Hence, it was decided that the DC flow resistance for each sample could be 
obtained from a least squares fit of the low frequency portion of the flow 


resistance results, as shown in Figures 23, 25, and 27, for that sample. 


The results are summarized in Table 5.3. 


Table 5.3 DC Flow Resistance 


‘aioe (cto arena 
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Therefore, theoretical values of the complex effective mass 
density can be calculated using Equation (5.18) for each of the porous 
solid samples. However, the quantities of interest are the effective fluid 


mass density and the flow resistance. 


Ze Effective Fluid Mass Density 
The (real) effective fluid mass density, Pog (@) , is found from 


Boy (Equation 5.18), by 





R fs 
Pog (@) = Rel Big] = py + a Im[F(x)] . (5.23) 


A plot of the effective mass density of the air contained in the 150 micron 
sample and its theoretical effective mass density calculated using Equation 
(5.23) is shown in Figure 29. Similar plots for the 300 micron and 500 
micron samples are shown in Figures 30 and 31, respectively. 

Although the experimental results for the 150 and 500 micron 
samples do not entirely agree with their theoretical results, there is good 
agreement between the experimental and theoretical results for the 300 


micron sample. 
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Figure 29 Effective mass density - 150 micron sampke. 
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Figure 30 Effective mass density - 300 micron sample. 


56 


density (kg /m-3) 
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Figure 31 Effective mass density - 500 micron sample. 


=, Flow Resistance 
The frequency dependent flow resistance, Riow(®) , is found 


from Puy (Equation 5.18), by 


Rygy(@) = - @Im[ 6,7] = Roc Re[F(x)] . (5.24) 


A plot of the flow resistance of the air contained in the 150 micron sample 
and its theoretical flow resistance calculated using Equation (5.24) is shown 
in Figure 32. Similar plots for the 300 micron and 500 micron samples are 


shown in Figures 33 and 34, respectively. 
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Figure 32 Flow resistance - 150 micron sample. 
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Figure 33 Flow resistance - 300 micron sampte. 
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Figure 34 Flow resistance - 500 micron sampke. 


Although the experimental results do not entirely agree with the 
theoretical results, especially at the high frequency limit, they do show a 
similar trend. Baker points out that the high frequency limit of flow 
resistance should be proportional to w! and that the low frequency limit of 
flow resistance should be proportional to wo, This is clearly the case, as 
shown in Figure 35. The log scale was chosen since it makes the 


calculated limits appear as straight lines and thus makes the comparison 


easier to visualize. 


59 


x — 150 micron sample doto 
+ — 300 micron sample data 
o — 800 micron sample dota 


+ ¢ t+ eeete 


Rftow (kg/s /m-3) 





103 
101 102 105 
frequency (Hz) 


Figure 35 Flow resistance - 150, 300, and 500 micron samples. 
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Vi. CONCLUSIONS AND RECOMMENDATIONS 


An experimental apparatus was designed and built to obtain the 
effective mass density and flow resistance of air in a porous solid from 
acoustic pressure and acceleration measurements. 

Initial experiments were conducted to test the apparatus and 
equipment setup on a known quantity, namely air in an empty sample tube. 
The objective of these experiments was substantially achieved after minor 
modifications to the apparatus and a revision to the experimental procedure 
were made. 

Porous solid sample experiments were conducted using three samples 
with mean pore diameters of 150, 300, and 500 microns. The results of 
these experiments were encouraging, especially for the 300 and 500 micron 
samples. However, the flow resistance results had a much closer agreement 
with theory than the effective mass density results. All results were 
repeatable and consistent from sample to sample. 

Future research should focus on the 300 and 500 micron samples 
since they produced the best results. Areas of interest include the 
possible transition between Poiseuille flow and boundary layer flow at 300 
Hz vice the predicted frequencies, and the corrections required to make 


the theory fit the data. 
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APPENDIX A PROGRAM FOR INITIAL EMPTY TUBE CALCULATIONS 


initial empty tube calculations 


gravity, m/s/s 
g = 9.8; 


speed of sound in air, m/s 
c = 343; 


frequency 
freq = [50:10:490,500:50:1000]; 
omega = 2*pix*freq; 


wave number 
k = omega / c; 


accelerometer sensitivities, mV/g 
lasen - 4.75; 
rasen = 4.5 6.: 


pressure sensor sensitivities, mV/Pa 


right 
rpsen = 10°((16.13 - 0.0015 * freq) / 20); 
rpsen = - rpsen; 
left 
Ipsen = 10°((15.53 - 0.00035 * freq) / 20); 
lpsen = - l]psen; 


Note: negative sensitivities due to reverse polarity of 
both accelerometers and pressure sensors 


length of sample tube, m 
] = 2102; 


62 


dP gP oP 


oP dP dP 


load emty.dat 


lpv = emty(1:56,1) + j*emty(1:56,2); 

rpv = emty(57:112,1) + j*emty(57:112,2); 
lav = emty(113:168,1) + j*emty(113:168,2); 
rav = emty(169:224,1) + j*emty(169:224,2); 


pressure 
rp rpv ./ rpsen.'; 


lp lpv ./ lpsen.'; 
acceleration 

la = lav / (lasen/g); 

ra = rav / (rasen/g); 

ra = - Ya, 


Note: the accelerometers mounted to the surface 
of the speaker cones were calibrated such that 
the direction of positive acceleration is outward 


acceleration correction factor 
acf = (8.3 / 8.9)°2; 


finite wavelength correction factor 
fwef = tan(k*1/2) ./ (k*1/2); 


= rp - Ip; 
grad_p = delta_p / 1; 


accel = (la + ra) / 2; 

accel = accel * acf; 

result = - (grad_p ./ accel) ./ fwef'; 
density = real(result); 

Rflow = - omega’ .* imag(result); 
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de 


dP JP 


APPENDIX B PROGRAM FOR REVISED EMPTY TUBE CALCULATIONS 


revised empty tube calculations 


anew SS sw ewes ew ewes ee e—e= wwesss—_ ase Ss were Berge eees sw owes ww ee ST ewew =|] == = = == 


een eae awe ee ee eae ee ee CellC ETUC reThllUCEUEUlUCUEUCCUOEOUCCECU LUC EEOC ll Sel eee ese ee ee eet Cee teeta ee ee eee eee ee eee Oe eee eee eee eee 


gravity, m/s/s 
g = 9.8; 


speed of sound in air, m/s 
© = 13543- 


frequency 
freq = [50:10:490,500:50:1000]; 
omega = 2*pixfreq; 


wave number 
k = omega / c; 


accelerometer sensitivity, mV/g 
rasen = - 4.56; 


Note: negative sensitivity due to reverse polarity 
of accelerometer 


pressure sensor sensitivity, mV/Pa 
lpsen = 0.3017; 


am ewe eweee@ es se ae es eS a es es Se SS es eee es eee ee ee ee ee Cem ch ell emcee ele oe oe ee 2 


el 


effective length of sample tube, m 
1 = .092; 


ante es SF ee SF aes SF aT ea ae ee eee Be eae ese ee BP ee ee ee Se SET ee eSB SE ee ee eee eee ee Ee eel Uc Ee UlUeeerml ee eee ee ee Se 


load empty.dat 
empty(1:56,1) + jxempty(1:56,2); 
empty(113:168,1) + j*empty(113:168,2); 


— 
im? 
<¢ 
i ott 


64 


dP dP oP 


dP oP de 


pressure 
lp = lpv / lpsen; 


acceleration 
ra = rav / (rasen/g); 
ra = -ra; 
Note: the accelerometer mounted to the surface 


of the speaker cone was calibrated such that 
the direction of positive acceleration is outward 


acceleration correction factor 
acf = (8.3 / 8.9)°2; 


finite wavelength correction factor 
fwef = sin(k*1) ./ (k*1); 


delta_p = - lp; 

grad_p = delta_p / 1; 

accel = ra; 

accel = accel * acf; 

result = - (grad_p ./ accel) ./ fwcef'; 
density real(result); 


a 
rh 
— 
Oo 
= 
i“ ott 


- omega’ .* imag(result); 
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oP 


oP oP 


APPENDIX CC PROGRAM FOR POROUS SOLID SAMPLE CALCULATIONS 


porous solid sample calculations 


=_—aeanmnemeeeeeeeeqmrmwemenm =e =e we sw we eee amewewe=a= = we ween ww wesw aweewmPemeeweewee = awe a= a= == a= 


neae7egawe & es see eaeaq@wae ae Se see SS Set eee s&s se weweegee@w@e@ge @&ew@g @& @fe@weeiea «ea @ == ==» 


gravity, m/s/s 
g = 9.8; 


frequency 
freq = [50:10:490,500:50:1000]; 
omega = 2*pixfreq; 


accelerometer sensitivity, mV/g 
rasen = - 4.56; 


Note: negative sensitivity due to reverse polarity 
of accelerometer 
pressure sensor sensitivities, mV/Pa 
lpsen = 0.3017; 
cpsen = 0.2775; 


eaag@gre= @& ss eawwqwqmwqqwFw@e°u°wei@qewe@wees ase Se ee Se S&S oe eweeeeeees we wedge é@ée®@=@ée «= «¢ 


load samplel5.dat 


lIpv = samplel15(1:56,1) + j*samplel5(1:56,2); 
cpv = samplel5(57:112,1) + j*samplel5(57:112,2); 
rav = sample15(113:168,1) + j*samplel5(113:168,2); 


ae @& ewes we we wwe ws we we weseeweeewmwee Ge ew es ewe Bees Be eweeesewr TW Ss SS SS SS SB SP SP SB | == = 


EE~ ee ee ew @weemeeaee=s Ss we we we wewmweewerweeenss SB ewrenwreenwe ss BB we we BWP swswewe  s—F ss BP BP se S| — 


effective length of the sample, m 
1 = .092; 


porosity of sample 
P = 0.37; 
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=_—ea~=s @ ase ae ae ewe ee ewe ewe awe ewe wee we Se ee ee ee ee ee ee ee es eee see eee ee 


pressure 
lp = lpv / lpsen; 
cp cpv / cpsen; 


acceleration 
ra = rav / (rasen/g); 
ra = - ra; 


Note: the accelerometer mounted to the surface 


of the speaker cone was calibrated such that 
the direction of positive acceleration is outward 


acceleration correction factor 
acf = (1/P) * (8.3 / 8.9)°2; 


finite wavelength correction factor 


kl =e acos( Ip ./ (2 * cp) ); 
fwef = sin(kl) ./ (kl); 
final calculations 

delta_p = - lp; 

grad_p = delta_p / 1; 

accel = ra; 

accel = accel * acf; 

result = - (grad_p ./ accel) ./ fwef; 
density = real(result); 
Rflow = - omega’ .* imag(result); 
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